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Abstract 

We have previously studied properties of a one-dimensional potential with N equally 
spaced identical barriers in a (fixed) finite interval for both finite and infinite N. It 
was observed that scattering and spectral properties depend sensitively on the ratio c 
of spacing to width of the barriers (even in the lirait N — > oo). We compute here the 
specific heat of an ensemble of such systems and show that there is critical dependence on 
this parameter, as well as on the temperature, strongly suggestive of phase transitions. 

PACS number(s): 05.70.Fh, 02.60.Cb, 03.65. Ge 



1 Introduction 

We have studied the one-dimensional locally periòdic multibarrier potential of finite range for 
both finite and infinite number N of barriers. We found |I| that there is a critical dependence 
of the transmission coefficient, the cross section and the distribution of poles of the S'-matrix 
on the ratio c of the total interval between the barriers to their total width. Under certain 
conditions this model was found to contain signatures of chaotic behaviour fl3|, such as ràpid 
spread of wave packets and Wigner type spectral characteristics. 
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We discuss in this work thermodynamic properties || such as the specific heat and the 
entropy of an ensemble of such systems and show that they also depend critically upon this 
ratio in addition to their dependence upon the temperature. We show that when the number 
of barriers is not large there are no low-lying energy eigenvalues for small vàlues of the ratio 
c. These vàlues of c depend upon the value of the number iV of barriers in such a manner 
that the ranges of c in which this part of the energy spectrum is missing increase at first 
with N and then decrease as N continue to grow until they disappear entirely for large 
enough N. Above the upper boundaries of these ranges are the points where the energy 
spectrum contains eigenvalue distribution over all energies. A similar phenomenon occurs 
in the work of Fendley and Tchernyshyov || on one-dimensional phase transitions, where 
the seriously disordered phase behaves as if it were at very high temperature, where the 
effect of interactions is relatively small, i.e, as an almost independent particle model. This 
apparently accounts, even for neighbouring vàlues of c within these ranges, for large changes 
in the average energy and all the other statistical mechanics properties, such as specific heat 
and entropy, derived from it. For example, the corresponding curves of the specific heat, for 
these neighbouring vàlues of c, as functions of the temperature, differ markedly from each 
other, as will be shown in the following, when plotting these curves in a single figure. This 
sensitivity of the specific heat is largest for intermediate vàlues of N and c in which the 
curves of the specific heats as functions of the temperature differ to the extent that some 
of them exhibit double peak phase transitions while other curves, for neighbouring vàlues 
of c, resemble the conventional Debye curve H] which characterizes the sòlid crystal. For 
infinite N, as will be shown, the double peak phase phenomenon is retained even for large 
vàlues of c. Also, we find, for all vàlues of iV and c, indication of phase transitions in specific 
heat at small vàlues of the temperature T. We note that for any c, which is greater than 
some value which depends upon N, the system behaves for all N, except for a spike form 
at small T, in a manner similar to that of the sòlid crystals in which the constituent àtoms 
are widely separated, and thus, as we have remarked, are characterized by large vàlues of 
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the ratio c. This is seen on the corresponding graphs of the specific heats as functions of 
the temperature which have the same form as that of Debye (we do not imply here that our 
system is analogous to a system of lattice vibrations, but only to the behaviour of the result 
due to the apparent presence of uncorrelated modes). 

In Section 2 we study the properties of a one-dimensional N potential barrier system 
when N is a finite number. In Section 3 we discuss the limit of iV — > oo (in the same fixed 
interval) . For both cases we find abrupt and large changes in the vàlues of specific heat and 
entropy for small vàlues of the temperature T. These large changes suggest the existence of 
phase transitions. In particular, for intermediate N and c, we find double peaks [|, [7|, j|, ||] 
in the specific heat curves. Double peak phase transitions appear, especially, for infinite N 
in which it remains effective even for large vàlues of c. Tanaka et al |7] have found double 
peak structures in anti-ferromagnetic materials, corresponding to magnètic phases, where 
the external magnètic field seems to play a role somewhat analogous to the parameter c 
in our study. Leung and Neda || and also Kim et al || have found double peaks in the 
response curves, apparently associated with dynamically induced phase transitions. Ko and 
Asakawa |§ have also found a double peak structure in their calculations of the phases of 
a quark-gluon plasma, where one may think of a large number of interactions in a bounded 
region. We do not imply that our model contains an analog of their interacting systems, but 
suggest that some of the mathematical properties may be common. 



2 The one-dimensional N potential barrier system 

We consider a finite N barrier system where all these barriers have the same height v and are 
locally periòdic in the finite interval. This array is assumed to start at the point x = — ^ 
and ends at x = so that the total length of this system is L = a + b. Here a is the total 
width of all the N barriers (where i; ^ 0), and b is the total sum of all the intervals between 
neighbouring barriers (where v = 0). Thus, since we have N potential barriers the width of 
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each one is and the interval between each two neighbouring ones is j^zï- Denoting b = ac 
where c is a real number we can express a and b in terms of L and c as [[[] a = b = j^. 

Let us first consider the passage of a plane wave through this system, which has the form 



A x é* x + B x e 



-ikx 



X < 



a+b 



. Matching boundary conditions at the beginning and end 



of each barrier, we may construct a solution in terms of the transfer matrices ||, [T]] on 
the jth barrier. After the nth barrier we obtain, using the terminology in ||, the following 
transfer equation [ïj 



p{n)p{n-\) _ _ _p(2)p(l) 



where is a product of three two dimensional matrices (ï|. A 2n +i and B 2n+ i are the 
amplitudes of the transmitted and reflected parts respectively of the wave function at the 
nth potential barrier. A\ is the coefficient of the initial wave that approaches the potential 
barrier system, and B\ is the coefficient of the reflected wave from the first barrier. Eq ([]]) 
may be written as fl]] 
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where k is y^pr- Note that the transfer equation (§) is vàlid for both cases of e > v and 
v > e except for the two-dimensional matrix in T which, for the e > v case, assumes the 
form 
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and for the v > e case 



T n = cosh(— ) + - smh(— ) 

f] ag. 
?2i = -~smh( — ), T 22 
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The q, £ and 77 of Eq (0) are |I[ ç 
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í = f + |> = f - ^, and those of Eq 
= — z£ = — — . In the numerical part of 



this work we assign h — 1, and m = \. Note that the two matrices T from Eqs (|3|)-([§) do 
not depent on n. 

We find, now, the energy spectrum of this system. We use the 5-matrix and the periòdic 
boundary conditions at the two remotely placed sides of the system. That is, we assume 
that the wave function and its derivative at the far right end of the system, say at x = C 
where C is much larger than the size L = a + b of the system, are equal to the wave function 
and its derivative at the corresponding far left end of the system at x — — C, so we obtain 
A2N+ie lkC = Aie~ lkC , B 2 N+ie~ tkC = Bié kC . Thus, using the last relation and those between 
the components of the S and Q matrices @, where Q is given by Eq (|2|), 
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one may write the following dependence of the outgoing waves A 2 n+i and B\ upon the 
ingoing ones A\ and B 2 n+i 
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To obtain a non trivial solution for the vector 
equation 



A2N+1 



one must solve the following 
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The form of the equations (ff)-(|6|) are the same for both cases of e > v and v > e except 
for the two-dimensional matrix Q from Eq (||) in which the matrix T assumes the form 
(§) for the e > v case and the form @ for v > e. The energies, satisfying this relation, 
depend, as noted, on c and N and are obtained numerically. We show, now, that the high 
energy part (e » v) of the spectrum obtained from Eq (|6|) depends only on C and may be 
obtained analytically without using numerical methods. For that matter we note that when 
the energy e becomes very large we have k ~ q, and £ and 77 (see the inline equations after 
Eq (f|)) obtain the vàlues of 2 and respectively. In that case the components T í2 and T 2 i 
of the two dimensional matrix T from Eq @ become zero, and the diagonal elements Tu 

iqa iqa 

and T 2 2 become e~ and e~~ respectively. Thus, the two dimensional matrix Q from Eq 
(§) becomes much simplified and may be calculated analytically as follows (using q ~ k for 
the very high part of the energy spectrum) 
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for n = N, required for the application of (|6|). Thus we see that the two dimensional matrix 
Q becomes the unity matrix in the limit of large energies e. Substituting the resulting 
components of Q (Qn = Q22 = 1, Q12 = Q21 = 0) in Eq @ we obtain 



det 



3 2ifcC 



- 1 



e 2ikC - 1 



(8) 



The last equation is satisfied for all k — ^r, n = 0, 1, 2, 3 . . ., but since we are restricted 
here to the high energy part of the energy spectrum we refer only to the large vàlues of n. 

As remarked, the previous equations may be applied also for the e < v case except for the 
correct Q. Thus, the energy spectrum is composed of three parts: 1) the part for the e < v 
case. 2) the part that satisfies e no > e > v, where e no is some arbitrarily specified large energy 
(these two parts are obtained numerically by solving Eq @ for both cases of e > v and e < v 
case with the appropriate Q), and 3) all the high energies that are larger than e no and are 
obtained analytically using Eqs (0), @ and the relation k = n = n ,n + l,no + 2 . . ., 
where n is some specified large integer that corresponds to the energy e no . Now, since the 
relation between the energy e and k is k = y/e (as remarked, we assign U — 1 and m — h), the 

high part of the energy spectrum is given by e^gh = (77 ) 2 , where n — uq, no + 1, no + 2, 

For the numerical part of the calculations we assign the following vàlues for the relevant 
parameters: v = 60, C = 90, L = 20, and e no = 1080. Thus, the e < v part of the spectrum 
is < e < 60. The second part is 60 < e < 1080, and the third part is e > 1080. The 
n that corresponds to e no is n = ^yfe~^ ~ 941. The assigned vàlues for C and L yield a 
ratio of f = 4.5 which ensures that the points x = ±C are far enough from the region of 
the potential and so we may assume (see the discussion after Eq @) that the vàlues of the 
wave function and its derivative at x — +C are equal to these at x — —C. These equalities 
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are needed for obtaining Eqs (|5|)-(|6|). Thus, the results of the numerical simulations depend 
upon the ratio ^ and not upon the specific vàlues of C and L. The somewhat higher value 
of the potential v = 60 (the value, conventionally chosen for numerical simulations (see |T3|| ), 
is from the range (2 < v < 8)) is especially chosen for the e < v part of the simulations 
in order to accumulate enough data and thus to obtain a better statistics. Note that the 
simulated e no > e > v part, where e no = 1080, yields a large amount of data so that in order 
not to remain with a comparatively small amount for the e < v part we have chosen, as 
remarked, a somewhat large value of v. The value of e no = 1080 is chosen as a limit value 
beyond which the corresponding terms of the sums Ee(c,JV)e-^ and E e _/3e(c,JV) in the 
following Eq (||) may be approximated by their simplified analogs obtained using Eq (§). 
Thus, the results are not sensitive to these specific vàlues of C, L, v, and e„ . 

We want to obtain a formula for the average energy from which we may derive most of 
the statistical mechanics parameters mentioned above. This average energy depends upon 
c, C, N and the temperature T (the dependence upon the temperature is through 8 — -4^, 
where k is the Boltzman constant which is assigned, in our numerical work, the value of 
unity) and is given by 

<e> (3 c N) — ±±_X2_Z_ +^94il ç) e 

where the first sum in the numerator and denominator contains the contributions from all 
e(c, N) < 1080. For higher vàlues of e the expressions are simpler and we take this into 
account in the second sum over all integers 941 < n. The first sum in the numerator and 
denominator of Eq (|]) includes the energies from the e < v and the v < e < 1080 parts of 
the spectrum . These parts are obtained numerically from Eq @ in which we substitute for 
the components of Q from Eq @, using the T's of Eq (||) for the e > v case, and those of Eq 
(|) for e < v . From Eq (||) we obtain an average energy for each specific triplet of vàlues for 
c, N, and 3, and from this average energy we may derive the quantities of statistical physics. 
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We note that the sums over the energies from the range 0.1 < e < 1080 depend upon the 
parameters N and c whereas the sums over the higher energies do not depend upon them 
(see Eqs (|7]), (|8|)). The specific heat Ch is obtained as the derivative of the average energy 
from Eq @ with respect to the temperature T. 

It is found that for large vàlues of the temperature T the curves of the specific heats Ch, 
for all vàlues of c and N, tend to the constant value Ch = 0.55. That is, for these T's the 
curves of Ch become as expected a constant curve as for Dulong-Petit |§]. Also, for small 
T's the curves of Ch, for all c and N, rise rapidly to their maximum vàlues Ch max from which 
they descend either to the asymptotic value of 0.55 for large T as noted or to some minimum 
from which they rise again to a second peak that descends to the value of 0.55 for large T. 
We note that the Ch max s are points at which the derivative appears to be very large and 



are, therefore, suggestive of the existence of phase transition [Ï0|]. The vàlues of these Ch max , 
however, as well as the behavior of the specific heat for intermediate vàlues of T depend upon 
c and N. It has been found that for large c the curves of the specific heats, as functions of T, 
are of the Debye type ||], that is, the ràpid approach to maximum Ch max for small T and the 
immediate decrease to a constant value as T grows. For small c, however, the forms of the 
specific heat Ch for intermediate T depart markedely from that of Debye, and the range of 
T in which Ch is different depends upon N in such a manner that this range increases as N 
grows. For example, when N = 2 this range is 0.4 < T < 3.8, for n = 35 it is 0.4 < T < 100 
and for N = 250 this range increases to 5 < T < 2000. 

Figure 1 shows 38 curves of the specific heats, for N = 6, as functions of the temperature 
in the range 0.1 < T < 35. Each curve is for a different integral value of c in the range 
2,3,..., 39. The central and dense part of the figure, where a large part of the curves have 
the same form, are those graphs that have a large c and, therefore, may represent the sòlid 
crystals that are characterized by a periòdic structure in which neighbouring occupied sites 
are widely separated. Indeed, these curves resemble, except for the sharp peaks, that of 
Debye which represents well the sòlid crystal. The other curves that differ from the central 
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ones, and that generally have large vàlues for the specific heats at small T, are those that 
have smaller c and, therefore, do not have the behaviour of the specific heat curves of a sòlid 
crystal. These curves may represent some "soft" substance (T^| in which the constituent 
àtoms or molècules are closer to one another than in the sòlid crystal. These substances 
have Einstein freqüències [|, [Ï2|] smaller than those of the sòlid crystals by a factor of 10 



to 50 |12| and, therefore, are characterized by higher vàlues of the specific heats. Also, 
the remarked sharp peaks have a very large value for the derivative and this suggests an 
existence of phase transitions. We have, especially, studied the immediate neighbourhood 
of the region containing ràpid variation in specific heat as a function of temperature and 
the parameter c (i.e, the neighbourhood of a phase transition). We see that there is a 
very strong and critical dependence on c. This result illustrates the physical mechanism 
for the formation of such ràpid transitions, a resonance-like phenomena controlled by the 
geometry of the barrier system. This is demonstrated in Figure 2 which shows 15 curves of 
the specific heat, for = 6, as functions of the temperature T and for the following vàlues of 
c = 0.3, 0.4, 0.5, . . . 1.7. The central dense part of the figure, which is composed of 9 curves, is 
drawn for the larger vàlues of c whereas the remaining 6 curves are for the smaller c's. Note 
that although the difference in c for any two neighbouring curves is only 0.1 nevertheless the 
upper curves in the figure differ significantly in appearance from each other and from those 
of the central part. 

Fendley and Tchernyshyov [||] have discussed one-dimensional phase transitions in sys- 
tems with infinite number of degrees of freedom per site. They argue that a singularity of 
the maximum eigenvalue in the "transfer matrix" (their model considers a set of systems 
with SU(N) type symmetry) causes a phase transition. The iV — > oo limit of || correspond- 
ing to an infinite number of degrees of freedom per site is replaced, in our case, by a large 
(infinite) number of barriers in a finite interval. The transfer matrix that we have introduced 
connecting neighbouring sites (barrier-gap structures) is independent of j3 (inverse temper- 
ature), but the eigenvalue of the total transfer matrix of the infinite system has a branch 
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point as we show in Section 3. This singularity can influence the behaviour of the partition 
function, resulting, as in ||], in the phase transitions that we observe in our numerical study. 
Figure 3 shows 39 different curves of the specific heat Ch as function of the temperature for 
N = 15. Each curve is for a different integral value of c from the range 2,3,4, . . .40. As 
in Figure 1 one can see a dense batch of similar curves in the central part of the figure and 
other 16 curves that are graphed one above the other in the upper part of it. The dense 
batch corresponds, as in Figure 1, to the larger vàlues of c and so may represent the sòlid 
materials that are characterized by a periòdic structure (see the discussion of Figure 1). The 
other 16 curves correspond to the smaller vàlues of c but compared to the former figure (for 
N = 6) one can see that these curves demonstrate the double peak appearance found in 
antiferromagnetic and superconducting materials ||. The first peak may be clearly seen 
at low temperatures at about T ps 0.5 and the second higher peak at about T ps 5. All 
the curves of Figure 3, the Debye-like singly peaked as well as the doubly peaked curves, 
merge together for large T into one batch that tends to the value of Ch(T » 1) ps 0.55. 
Note that Figure 1 also shows the same general form of a central dense batch of Debye-like 
curves and other different graphs in the upper part of the figure. These graphs, however, 
show no sign of double peak, even when finely graphed in the neighbourhood of the critical 
temperature. Thus, one may conclude that the double peak phenomenon is related to the 
number of barriers so that it is more apparent for the large number of them as seen from 
Figure 4 which shows 39 different curves of the specific heats, for n = 35, as functions of 
the temperature in the range 0.1 < T < 100. Each curve is for a different integral value 
of c in the range 2,3, .. . ,40. As in the former figures the similar Debye like curves in the 
central part of the figure are for large c's and the partition function, for these vàlues of 
the parameter c, therefore should contain some mathematical features in common with the 
partition function for the vibrational modes of sòlid crystals. The other curves are for small 
c's and they may correspond , as in the former figures, to "soft" substances. We note that 
seven of the curves have each a part below the Debye-like curves and a part above them 
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and so they resemble, in a more apparent manner than Figure 3, the remarked double peak 
phenomena |7|, ||. The second peak, in our case, is obtained at a comparatively large value 
of the temperature (T « 40) compared to the vàlues (T « 1) in || [7|, |P and to the value of 
T « 5 in Figure 3. Note also that the maximum vàlues obtained by the curves of this figure 
are unity, whereas, most of the curves in the former figures have màxima that exceed unity. 
All the curves of figure 4 show for small vàlues of T, as in the former figures, peaks that are 
suggestive of phase transitions. 

Thus, as remarked, the double peak appearance is more pronounced for the larger vàlues 
of N, but we note that as N becomes larger the second peak diminishes in height until it 
disappears entirely for large enough N and remains only the first peak for small T. But, 
as we see in Section 3, for the limit N — > oo the double peak phase transition is seen for a 
larger range of c than for finite N. 

We can find the corresponding critical exponent |Ï0| x associated with these phase transi- 
tions by noting, after studying and analysing the behaviour of these curves in the immediate 
neighbourhood of the critical temperature T c , that we can write an analytical approximate 



expression for the specific heat, in the neighbourhood of these points, as follows |10 



C h (e) = A + Bev, (10) 

where e = T ^ Tc . As seen, the first order derivative of this specific heat with respect to the 
temperature diverges at the point T = T C and so, the critical exponent x is obtained as (Ï(J 

x = 1+lim M^M = 1 + lim ^ = i, 

A e-o ln(t) e^o l n \e\ 2' y ' 

where Ch{t) is the derivative of Ch from Eq ( |Ï0| ) with respect to e and the unity value of 
the first term denotes the order (which is 1 here) of the derivative of Ch from Eq ( |Ï0| ) which 
diverges at T = T c , that is, the appropriate critical exponent is |. 

The conspicuous departure of the curve of the specific heat, for small c, from that of 
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Debye-like behahiour can be explained by noting that the total number of nondegenerate 
energies, as in [|], that satisfies Eq (|6|) varies for different vàlues of these c's. Moreover, 
we íind, numerically, for all íinite N and for small c, no energy from the lower part of the 
spectrum that satisfies Eq (|) for the e < v case. That is, solutions of Eq (|6|) for this case 
are found, for small c's, only from the part of the spectrum that is close to the value of v. 
Thus, when we sum upon all the allowed energies, in order to calculate the average energy 
and the specific heat, the summation does not include the lower part of the spectrum. For 
example, the number of nondegenerate allowed (energies) solutions of Eq (H), for N = 6, 
e < v and c = 1.5 are only 3, whereas they amount to 311 for c = 15. That is, for the larger 
vàlues of c we have a larger number of additional (that may be thousands for large potential 
v) allowed energies, and this yields entirely different vàlues for the average energy < e > and 
the specific heat Ch derived from it. 

The entropy S and the specific heat Ch are related by Q Ch = Thus, from the last 

discussion we infer that also the change of the entropy S with the temperature T jumps at 
the same vàlues of T in which the specific heat Ch jumps. That is, the change of the entropy 
with T has also phase transition. Moreover, the entropy S and the free energy F are related 
by the equation [|] S = — so that, differentiating both sides of the last relation with 
respect to the temperature T and using the relation between the specific heat Ch and the 
entropy S we obtain 

„ ^dS d 2 F , 

= &ï = ~ T im (12 > 

From the last equation we see that the second derivative of the free energy F with respect 
to the temperature T also changes steeply at the same vàlues of T in which Ch does so. 

When c becomes very large we have b » o, so that we may ignore a compared to b. 
Thus, writing the trigonometric functions of the components of T from Eq (||) as exponentials, 
substituting in Eq (||), and ignoring, as noted, a with respect to b we can see that the matrix 
Q from Eq @ becomes in the limit of a very large c the two dimensional unit matrix. In 
this case Eq @, from which the energy spectrum is obtained, becomes the same as Eq (H) 
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from which the energy spectrum has been obtained as e = k 2 = ( z ^-) 2 . But we note that 
whereas Eq @ was obtained for the case of high energies only for which the index n begins 
from a large value n (no = 941), here, in the limit of very large c, n assumes all integer 
vàlues of n = 0, 1, 2, . . .. In this case the average energy is 

y^n=oo/7m\2 p -/3(^) 2 

< e>Qu. W = Vn=~ l-mY (13) 

Z^n=0 e 

Note that < e >c h ^ (P) from the last equation does not depend upon the number of barriers 
N. The specific heat is 

g < e > Ch _ (g) 1 Eg°(f ) 4 e-^^) a S"(f) 2 e-^^) 2 E|g^)V^g 

hc - >0 ° <9T T 2 V^e-^) 2 V n=«D -/3(^)a V n=oo -/9(#)a 

(14) 

Plotting the curve of the specific heat from the last equation as a function of the temperature 
(not shown here) one can see that at small T varies rapidly from zero to 0.52 from 

which it descends sharply to its asymptotic value of 0.5. The curve is not differentiable at 
the point at which it assumes the value of 0.52 and so this point appears to be a phase 
transition one. 

3 The one-dimensional N potential barrier system for 

N ^ oc 

We discuss, now, the case where the number of barriers N tends to the limit N — ► oo. We 
may use for this case the equations ([]])- © derived for the finite N case in the previous 
section, so that taking the limit of a very large N one obtains from Eq @ for the right hand 
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side of the potential barrier system at the point x = ^ where ra = iV (ÏJ 
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(15) 



The middle expression was obtained by expanding in a Taylor series the cosine and sine 
functions and keeping only terms of the order and the last result by using the relation 
lim Tl _ +00 (l + = e c , where c is some constant. The <t 2 , and a 3 are the two dimensional 

o -?; 10 

The second exponent of the last result of 

Eq ( ]T5D may be expanded in a Taylor series, so that after collecting corresponding terms we 
obtain 

(16) 
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where / and d are defined as ffl 



6 V 
f = kb + aq-, d = aq- 



Thus, making use of the relation (/cr 3 + ida 2 ) 2 = f 2 — d 2 = 2 , and defining z = k(a + b) we 
obtain from Eqs (|T5|)-(|Ï6|) for the e > v case QTJ] 
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(17) 
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For the e < v case we use Eqs (M, and the corresponding quantities M /, d and 



f=kb-^L, d=^, ^ = (fa 3 -ida 2 ? = f*-d*, 



to obtain the following matrix equation equivalent to Eq (fÏ7|). 
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e- iz (cos ó + ihím) _ ie -izàtà^) 



(cos (f) — i f b ™^- > ) 



A 1 
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(18) 



We can, now, find the energy epectrum of the dense system in an equivalent way to the finite 
N case of the previous section. For both cases of e > v and e < v , we obtain equations 
similar to Eq (|^), but now the two dimensional matrices Q are those on the right hand 
side of Eqs (|T7|) , (Pj|) , where their four components are given explicitly. Note that the four 
components of the two dimensional matrix Q for finite N (see Eq (|])) can be obtained only 
numerically. Thus, using, for the e > v case, the explicit expression of Q from Eq ([Ï7|) we 
can write the analogous equation (for N — > 00) to Eq (|6|) as 



det 



Q22 

è 



Q2 

2ikCQ 21 £ 2ikC 



Q22 



QÏ2 



2e 2ikC | i | Q 12 Q 2 ie AlkC 
Q22 Q22 



(19) 



, II ^ W1 d 2 sm 2 (ó). , .. 9 . ,. f 2 sm 2 (ó). 2 f sm(ó) , . . . , ,, 
= cos(4/cC)(l + — ^) + cos(2z)(cos 2 (0) - — ^) + - - yr > (sin(2z) cos(0) - 

-sm(2kC + z)) -2cos(2kC + z) cos(0) + ï(mn(4fcC)(l + -^) + sin(2^)(cos 2 (0) - 

_ / 2 sm 2 (0) ) + 2/sm(0) (cos(2fcC< + ^ _ _ ^ + ^ cog ^ = Q 



For the e < v case we use the explicit expression of Q from Eq (|Ï8l) to obtain a similar 
equation to Eq ( |Ï9^ ) from which the energy spectrum for the e < v case may be obtained. 
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Defining the parameters k and r as 

one can obtain the eigenvalues of either Eq (|Ï7|) or (|Ï8|) in the form 

Ai, 2 = tcos(0 — k) ± ^/r 2 cos 2 (0 — k) — 1 (20) 

The eigenvalues of the e > t> case are obtained by substituting the correct (j) (see the inline 
equation after Eq flÏ6|) ) and those of the e < v case by substituting the corresponding quantity 
(see the displayed equation after Eq (|Ï7|)). From Eq (|20|) one can see that the derivatives of 
the eigenvalues Ai 2 may be singular at certain vàlues of <$> so that they fulfil the condition 
in for finding phase transition in a one dimensional system. This condition is necessary 
but may not be sufficient. Our numerical results suggest that such a transition occurs. 

The energy spectrum is composed from those energies that satisfy the real and imaginary 
parts of the last equation for the e > v case and the corresponding one for the e < v case. 
Thus, we may obtain the average energy < e >jv-kx> (c) for each value of c as in Eq (g) 
(without, of course the dependence upon N). From these average energies we obtain the 
corresponding specific heats Ch(c) n-^oo as functions of c and T. Note that although Eq (|Ï9|), 
from which we derive the average energy <e>7v_» 00 (c) and the specific heat Ch(c) n^ooi is 
obtained analytically compared to the corresponding Eq (||) for finite N, nevertheless, these 
Ch(c) n->oo have the same form as those of the finite iV (obtained by differentiating Eq (||) 
with respect to T) and are therefore difficult to study analytically. This is true, especially, 
for small T where the phase transitions are generally encountered. 

The dependence of Ch(c)N^oo upon the temperature T, as a function of c, is, for small 
T, different from the dependence discussed in the previous section for finite N. That is, it 
jumps up to its peak value from which it immediately jumps down to rise again to anther 



3 THE ONE-DIMENSIONAL N POTENTIAL BARRIER SYSTEM FOR N -> oo 18 



higher maximum. We note that, generally, for finite N there is only one peak maximum, 
and although in Figures 3 and 4 we see several curves that have double peaks, nevertheless, 
this is only for small c and that when c grows the curves become the same as that of Debye 
as seen in the dense central part of Figures 3-4 which are for large c. Compared to this the 
double peak appearance of the specific heat curves for infinite iV is retained even for large c 
as can be seen from Figure 5 which is drawn for c = 200. Figure 6 shows 30 different curves 
of the specific heats as function of the temperature. Each curve is for a different value of c 
from the set 0.3, 0.4, . . . , 3.2. Note the large difference in the heights of the two peaks of each 
curve, and that both are points where the first derivative with respect to the temperature 
T attains a very large value and so they appear to be phase transition points. A similar 
discussion to that of the finite N case (see also Eqs ([Ï0|) , (|ÍT|)) yields a critical exponent of 
\ for both peaks. That is, producing these curves in a fine grained manner in the immediate 
neighbourhood of the critical temperature T c we realize, as for finite iV (see the discussion 
before Eq (110)) that we may approximate analytically the form of C/j(c) jv-^oc by Eq ( pT0|) 
with a critical exponent of \ ~ 0.5. Indeed, comparing the forms of the curves for finite N 
to those of the infinite N in the neighbourhood of the critical temperature T c one does not 
find a large differece. Moreover, as c grows both kinds of curves show generally, except for 
the double peak appearance which is retained for the infinite iV even for large c, the same 
behaviour which characterizes the Debye-type curves. This may be seen from Figure 7 which 
shows 40 different curves of the specific heat as functions of the temperature T for integral 
vàlues of c from 2,3,4,... 41. Comparing these curves to those of Figure 6 which are drawn 
for smaller c one sees that the dense batch of similar curves in the central part of the figure, 
which characterizes the sòlid crystals (as we have encountered in the former figures), have 
appeared also for the infinite N case. This is because of the larger vàlues of c, for which 
Figure 7 is graphed, that enable, as for the finite iV case, a Debye-like forms for these curves. 
Note that this form is absent in Figure 6 (and also in Figure 2) because all the curves there 
are drawn for small c. The other curves in Figure 7 that are not part of the central dense 



4 CONCLUDING REMARKS 



19 



batch are for the smaller vàlues of c and, therefore, they are similar to those of Figure 6. 
As noted, the difference between the finite and infinite iV lies, especially, in the double peak 
phenomenon that is seen in the infinite N case even for large vàlues of c as can be seen from 
Figure 5. When Figure 7 is produced in a fine grained manner in the neighbourhood of the 
critical temperature T c one can see clearly the double peak for any curve as in Figure 5. As 
the temperature increases all the curves tend to the value of 0.55 as for the finite iV case. 
When c becomes very large the curves (not shown here) of the specific heats become similar 
to each other and to the Debye graph. 

We infer from the former results that also the first derivative of the entropy S and the 
second derivative of the free energy F, both with respect to the temperature, change in an 
abrupt manner at the same vàlues of T (see the analogous discussion at the previous section). 



4 Concluding Remarks 

We have shown that the one-dimensional multibarrier potential of finite range shows signs 
of phase transitions in specific heat for certain vàlues of the temperature T. These phase 
transitions depend upon the number of barriers N and the ratio c of the total spacing to their 
total width and are demonstrated for both cases of finite and infinite number N as shown in 
figures 1-7 and also for small and large vàlues of c. Moreover, it is seen from the curves of 
the specific heat as a function of the temperature for N = 35 and small c (see Figure 4) and 
for infinite N and a large range of c (see Figures 5-7) that the phase transitions appear in a 
double peak form. Double peaks have been seen in antiferromagnetic and superconducting 
materials and are apparently associated with dynamically induced phase transitions || [?]] 
and in the quark-gluon plasma |J. 

We note that we have found HJ that the one-dimensional multibarrier system discussed 
here demonstrates also, for large N, a unit value for the transmission probability and signs 
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of chaos which may be interpreted in terms of effective decoherence and the space analog [Ï6 



of the Zeno effect ||14|| in which a very large number of repetitions of the same experiment 
(interaction), in a finite total time, preserves the initial state of the system. It has also 



been shown flT|, [Tj| that a beam of light that passes through a large number of analyzers 
arrayed along a finite interval of a spatial axis, a configuration which is very similar to the 
one discussed here, remains after the passage with the same initial polarization and intensity 
it had before passing. This kind of preservation of the initial "state" by passing through 
a large number of physical apparatuses, each of them is supposed by itself to change the 
state of the passing system, has also been shown in the classical regime |ï7j where the initial 



density of classical partides passing through a one dimensional array of imperfect traps [18 



remains at the same value it had before the passage if the ratio of the total spacing to width 
(which corresponds to the ratio c here) increases. 

Thus, our finding here that when c grows the curves of the specific heat, as functions 
of the temperature T, become similar to the known graph of Debye [|J, indicates that the 
system makes a transition to the physical situation which behaves like the vibrational modes 
of a sòlid crystal. 

We have found that the critical exponents associated with these phase transitions have 
the value \. The other statistical parameters associated with the specific heat such as the 
entropy and the free energy also demonstrate, in the rate of their changes with respect to 
the temperature T, the same type of behaviour at the same vàlues of N, c, and T. 
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Figure 1: 38 different curves of the specific heat as a function of the temperature for N = 6 
and the following integral vàlues of c = 2, 3 . . . 39. The dense central part of the graph are 
the curves obtained for the larger vàlues of c and they all resemble, except for the sharp 
peaks at the left, the Debye's graph for sòlid crystals. The curves have very large vàlues 
for the derivative at the sharp peaks and so may be suggestive of phase transition at these 
points. The other curves are for small c and they represent the "soft" substances that have 
a comparable low Einstein freqüències and therefore a high vàlues of the specific heats. 
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Figure 2: 15 different curves of the specific heat as a function of the temperature for N = 6 
and the following vàlues of c = 0.3, 0.4, ... 1.7. The dense central part of the graph is 
composed of 9 similar curves obtained for the larger vàlues of c but since these vàlues are 
small compared to those of Figure 1 these curves are not similar to the Debye graph. The 
other 6 curves are for the smaller vàlues of c and they deviate significantly from each other and 
from those at the center. Note that a small difference of only 0.1 in c between neighbouring 
curves is capable of producing this large difference in appearance. Comparing this figure 
with the former one which has been drawn for the same N but for larger vàlues of c one sees 
that only 3 curves out of 38 in Figure 1 deviate markedly from the dense central part. 
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Figure 3: 39 different curves of the specific heat as a function of the temperature for N = 15 
and for the following integral vàlues of c = 2, 3 . . . 40. As in figure 1 the dense central part 
of the figure are the curves that represent the Debye curve, except for the sharp peaks at 
the left at which the specific heats have very large vàlues for their derivatives and so are 
suggestive of phase transition. Note that the curves that are not part of the central dense 
batch have double peaks where the first one may be seen at low T. 
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T 

Figure 4: 39 different curves of the specific heat as a function of the temperature for A^ = 35 
and for the following integral vàlues of c = 2, 3 . . . 40. As in the former figures the dense 
central part of the figure are the curves that represent the Debye curve, except for the sharp 
peaks at the left at which the specific heats have very large vàlues for the derivatives. Note 
that the double peak character of the curves that are separate from the central part is more 
pronounced than at the former figure. 
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Figure 5: Double peak of the specific heat as a function of the temperature for infinite N 
and for c = 200. Note the difference between the heights of the two màxima and that the 
double peak is retained even at this large value of c compared to the finite N case. 
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Figure 6: 30 different curves of the specific heat as a function of the temperature for infinite 
N ans for the following vàlues of c = 0.3, 0.4, . . . 3.2. All the curves tend to the value of 0.55 
for large T. The first peaks are not clearly shown for all the curves due to the large range of 
T over which they are drawn. For smaller ranges of T, as in Figure 5, all the first peaks are 
clearly shown (not here). Note that the second peak of each curve is much larger than the 
first and that all these second peaks are obtained for T > 5. The peaks of each curve have 
very much large vàlues for their derivative and so they are suggestive of phase transitions. 
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Figure 7: 40 different curves of the specific heat as a function of the temperature for infinite 
N ans for the following integral vàlues of c = 2,3, . . .41. Note the central dense batch of 
the similar figures known from finite N. It appears due to the large vàlues of c compared to 
these for which Figure 6 (and also Figure 2 for the finite N) was drawn which is the reason 
that the Debye-like form is absent in All the curves of Figure 6 (and Figure 2). Producing 
the curves of Figures 6 and 7 in the immediate neighbourhood of the crtical temperature T c 
demonstrate clearly the existence of their first peaks (not shown here). 



